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Abstract. Protocols used in quantum information and precision spectroscopy rely on 
efficient internal quantum state discrimination. With a single ion in a linear Paul trap, we 
implement a novel detection method which utilizes correlations between two detection events 
with an intermediate spin-flip. The technique is experimentally characterized to be more 
robust against fluctuations in detection laser power compared to conventionally implemented 
methods. Furthermore, systematic detection errors which limit the Rabi oscillation contrast in 
conventional methods are overcome. 
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1. Introduction 

A standard protocol for quantum mechanics experiments is comprised of a three-step 
sequence: Firstly, the quantum system is initialized to a known quantum state, for instance by 
optical pumping in the case of atoms or ions. Secondly, external control fields, e.g. sequences 
of laser light pulses or radio-frequency radiation, are applied to implement the Hamiltonian 
of interest. Finally, the resulting state is read out by a projective measurement. In 
the case of single or few-atom quantum systems, repeating the experiment many times 
is necessary in order to accumulate statistics, such that the average measurement result 
represents the quantum state after the evolution according to the Hamiltonian, projected onto 
the measurement basis. Thus, the fidelity of state detection plays an essential role in the 
understanding of the underlying dynamics and, in addition to other constraints given by the 
experimental setup, imposes a limitation to the quality of the experimental findings. 

Trapped ions are an ideal system to investigate quantum mechanical effects owing to the 
almost perfect control of the ion's internal and motional state preparation, manipulation and 
read-out [1J. This is demonstrated by both relevant theoretical and experimental progress, 
providing new insights in many fields of physics, to name a few: quantum computing [0, 
optical clocks 01 HI, precision spectroscopy flSl [6l [Vj and quantum simulations lf8ll9l [T0l[Tni . 
All underlying experimental schemes depend on the high fidelity of both the applied algorithm 
and the state detection of the ion. 

The basic principle of state discrimination in ions relies on electron-shelving [Tl2l . Two 
different energy states (qubit) are distinguished by their state-dependent fluorescence via 
coupling to a third level. In its simplest form, the number of collected photons during a 
single detection cycle determines whether the ion is assigned to a so-called bright (dark) 
state depending on this number being higher (lower) than a chosen threshold. This threshold 
detection technique has been successfully implemented in many experiments yielding an 
almost unity read-out fidelity for qubits stored in two optically separated meta-stable states 
( 40 Ca+ [13], lu Cd+ 031, 88 Sr+ 031). Detection fidelities can be further improved by 
taking into account photon arrival times, as opposed to only observing the integrated signal. 
This type of Bayesian inference or maximum likelihood detection has been successfully 
demonstrated in an optical qubit with 40 Ca+ yielding fidelities of 99.991(1)% OH9- The 
detection fidelity is ultimately limited by the finite overlap between the bright and dark state 
photon distributions, resulting in an increased and asymmetric (bias) error in the assignment of 
qubit states. Long detection times reduce this error, but are ultimately limited in real physical 
systems by the finite lifetime of the qubit states due to spontaneous or induced transitions 
to other states. This effect is particularly prominent for hyperfine qubits, where the detection 
fidelity suffers from depumping of one of the qubit states ( 9 Be+ 031 ED, 171 Yb+ 0H, 25 Mg+ 
11201 ). effectively limiting the detection time and thus the number of observed photons. This 
can be overcome by implementing a quantum non-demolition measurement in which the qubit 
state to be detected is repeatedly transferred to an auxiliary ion, where detection is performed. 
A Bayesian state inference of the time series of photon detection events of such an experiment, 
implemented with an optical qubit in Al + , and read out by a hyperfine qubit in Be + , yielded 
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a detection fidelity of 99.94% II2TTI . However, achieving such a high fidelity requires either an 
auxiliary ion to implement a quantum non-demolition measurement [21], or an optical qubit 
with a long excited state lifetime. Here, we demonstrate a novel state detection technique 
with improved state discrimination by combining two detection events of either threshold 
or Bayesian detection types with an intermediate well-controlled state inversion spin-flip. 
Observation of the correlated detection outcome acts as a post- selective statistical filter, 
significantly reducing the state assignment errors due to overlapping photon distributions. 
Owing to the symmetry of the ^-detection scheme, the detection errors for the two qubit 
states are at the same time equalized, effectively reducing the bias error. Another implication 
of this effect is the robustness of the scheme against fluctuations of detection parameters, 
making it particularly well-suited for systems in which only very few photons are detected. 
Consequently, the ^-detection scheme will always improve the detection fidelity of threshold 
and Bayesian detection schemes if the state inversion can be performed with sufficient fidelity. 
The scheme is also applicable to multi-particle systems, as long as individual atoms can be 
resolved. Furthermore, we imagine similar schemes to be applied to neutral atoms with non- 
destructive detection methods where the particle is not lost after a single detection cycle, such 
as measuring the state-dependent reflection and transmission signals of molecules or atoms in 
cavities 112211231 

2. Experimental Apparatus 

The experiments discussed here have been carried out in a setup described in detail in 
ll20l . Briefly, we trap a single 25 Mg + ion in a linear Paul trap at trapping frequencies 
G) m ~ lit x 2.2 MHz in the axial and CD, ~ lit x 4.5 MHz in both radial directions. The 
qubit is encoded into two long-lived hyperfine ground states of the Sj/2 manifold, namely 
\F, mp) = |3, 3} = 1 1) and |2, 2) = 1 1), representing the bright and the dark state, respectively. 
For all measurements, a magnetic field of typically B ~ 0.6 mT lifts the degeneracy between 
the magnetic sub-levels. A frequency-quadrupled solid-state fiber laser system provides light 
at 280 nm for laser cooling and detection using the cycling transition between the | \) and the 
1 4, 4) = \e) state of the P 3 / 2 manifold with alinewidth of Y ~ 2k x 40 MHz. The detection and 
Doppler cooling lasers are derived from the same source, their only difference being that the 
detection laser operates on resonance with the cycling transition, whereas the Doppler cooling 
laser is detuned by half the natural linewidth for optimal cooling performance. 

The resulting resonance fluorescence scattered from the ion is collected by both an 
objective and a parabolic mirror and subsequently focused onto a photo-multiplier tube for 
quantitative readout on a /is time scale. The application of radio-frequency at 1.789 GHz, 
which corresponds to the hyperfine splitting in 25 Mg + , provides coherent coupling between 
the qubit states. The radio-frequency is supplied by a quarter-wave antenna at a distance of 
approximately 12 cm from the center of the ion trap. An overview of all relevant levels and 
transitions is shown in Figure [H(a). 
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Figure 1. (a) Magnesium level scheme. Only relevant levels and transitions are shown. 
Doppler cooling (DC) is applied on the cycling transition between the S i /2 and P 3 / 2 manifolds. 
An additional o-repumper improves state initialization. The qubit levels | 4-) and | f) are 
coherently coupled by applying radio-frequency pulses (RF) at 1 .789 GHz. (b) Experimental 
protocol. Three Doppler cooling and repumping pulses initialize the ion to the | \) state. 
After that, a radio-frequency spectroscopy pulse is applied. The sequence is concluded with a 
resonant detection pulse (same as the DC beam). To obtain sufficient statistics, the sequence 
is typically repeated for 300-1000 times. 



3. Experimental Protocol and Quantum State Detection 

After the ion is loaded into the trap, the experimental sequence comprises three steps (see 
Figure [lib)) which are repeatedly applied to the ion: 

1. Preparation and Cooling. The ion is initialized by optical pumping to | \) by 
means of a Doppler cooling laser which is applied for 600 /is. An additional repumping laser 
guarantees that no population remains in the upper manifold {Si/ 2 F = 2). After Doppler 
cooling, the temperature of the ion is close to the Doppler limit of ~ 1 mK 11201 . 

2. Radio-Frequency Spectroscopy. A radio-frequency spectroscopy pulse is applied 
for a time T spec to drive magnetic field induced Rabi oscillations between | \) and | f) . Due to 
the small Lamb-Dicke parameter of the radio-frequency transition r\ ~ 10 7 , the transition is 
independent of the motional state of the ion, which makes Doppler cooling sufficient for these 
experiments. 

3. State Detection. The quantum state of the ion is read-out using different protocols, 
involving at least one detection laser pulse on resonance with the | ^) •<->• \e) transition. At the 
end of each experimental cycle the ion is projected into either the | \) or the | f) state. The 
ion's state is inferred by averaging over the results of many identical experimental cycles. 

If the ion is in the | \) state, application of the detection laser at saturation intensity 
will let it reside in the cycling transition and scatter photons with a rate of approximately 
r/4 ~ 6.425 x 10 7 s _1 , where T corresponds to the natural linewidth of the magnesium 
transition ll24l . On the other hand, if the ion is in the | f) state , the detection laser is detuned 
by the ground state splitting of A HF 2k x 1.789 GHz and photon scattering is suppressed by 

r 2 /2A 2 . 
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The collected number of photons are random variables. For the bright state, the random 
variable follows a Poissonian distribution [|25l . The probability to measure k photons is given 
by 

where | = R ■ x denotes the average photon number, R is the scattering rate and x is the 
detection time. 



3.1. State Detection Errors 

There are several sources of systematic errors which yield false state information. They 
impose limits on the detection time and the read-out fidelity. In the following, these errors 
are briefly described. Afterwards, different detection methods which help overcome these 
limitations are discussed. A detailed theoretical discussion on detection errors of the threshold 
and Bayesian methods is found in 11261 . 



3.1.1. Depumping and Decay of the Upper Qubit State. In the case of a qubit encoded into 
two hyperfine ground states, as for 25 Mg + , the finite detuning of the detection laser pulse with 
respect to the | f) state results in off -resonant excitation of the ion with a rate of approximately 

^ ~ 2k x 2.7kHz (j ~l) , (2) 

where so = I/I s is the saturation parameter, / the laser intensity and I s the saturation intensity 
of the Doppler cooling transition in 25 Mg + . Consequently, for long detection times, the ion is 
off -resonantly pumped into the | \) state, where it enters the cycling transition. Once there, a 
significantly larger amount of photons is scattered, which leads to a wrong assignment of the 
ion to the bright state, although it was initially dark. This results in a bias error of the state 
assignment lfT3l [16112111 . thus imposing a limit on the detection time. The effect of depumping 
on the number of scattered photons is depicted in Figure[2](a). The average number of detected 
photons for the bright state is 

|b(t)=/? b -t , (3) 

where R B is the scattering rate of the bright state. For the dark state, the depumping 
effect is taken into account in the following way: If the ion remains in the dark state 
during the detection interval, only background counts % x = R x ■ x are measured. Here, Roo 
corresponds to the background scattering rate and x to the detection time. In case the ion is 
depumped at a time t during the detection interval, the number of photons depends on the time 
(|; — Root + R B (x — t) see 11271 and Section I3T21) . Both processes are weighted exponentially by 
the function w(t) = exp(—t/T)/T, where T denotes the 1/e-decay time. The average number 
of detected photons for the dark state then reads 

| d (t)= j*dtw{t)l t + J dtw{t)loo = R*-x-(R B -Roo)T(\-e- xlT y (4) 
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The fit of the measured data yields a decay time of T ~ 61 ± 4 /is (see Figure [2a)), which is 
consistent with the expected depumping rate of Eq.[2] 

In principle, the bright state can also be depumped to the dark state if the polarisation 
of the detection light is not pure. However, this effect can be easily eliminated and will not 
be considered here. By comparing the initial bright state preparation in the optical pumping 
cycle with and without the c-repumper laser, we typically observe an estimated infidelity of 
1-2 % in the polarization of the detection light. 

One possibility to circumvent the depumping process is to map the hyperfine qubit in a 
state- selective way onto an optical qubit encoded in a ground and a meta-stable excited optical 
state. Since the magnitude of the detuning of the detection laser is in the optical regime, off- 
resonant depumping can be neglected. However, the finite lifetime of the excited meta-stable 
state also imposes a limitation on the detection time since the ion eventually decays back to 
the ground state by spontaneously emitting a photon. In addition to that, the finite fidelity of 
the mapping protocol has to be taken into account [fT3l . 




Detection Time t (us) No. of Photons 

Figure 2. (a) Depumping of dark state. The measured average number of collected photons 
as a function of the detection time is shown. The photon number increases linearly, if the ion 
is in the bright state (red circles). The solid red line corresponds to a fit to Eq.|3] For an ion 
in the dark state [black circles), the count rate reflects only the background photon counts for 
short times. The detection laser eventually depumps the ion to the bright state, thus the count 
rate asymptotically approaches the rate of the bright state for long detection times. The solid 
black line corresponds to a fit to Eq.|4]and yields a decay time of T « 61 ±4/is. The dashed 
green line reflects the asymptotic behaviour of the dark fluorescence for large detection times, 
(b) Measured normalized example photon histograms of the bright and the dark state for a 
detection time of 10 /is for a sing le 25 Mg+ ion. The measured scattering rates for the bright 
(dark) state are R B ~ 249 • 10V (/?,*, ~ 11 ■ lO^" 1 ). 



3.1.2. Overlapping distributions. In Figure [2](b), the measured probability distributions of 
the collected fluorescence photons for both bright and dark states of the 25 Mg + ion in 10 /is 
are shown. As a consequence of the limited detection time and the depumping events, the 
distributions exhibit a finite overlap. If the number of detected photons lies in the overlapping 
region, it cannot be unambiguously decided to which state the number of collected photons 
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should be assigned, leading to false state assignments with a bias that depends on the fraction 
of the photon histograms below (above) the threshold for the bright (dark) state. 



3.1.3. Fluctuations of experimental parameters. Temporal fluctuations in the atomic 
scattering rate induced by intensity or frequency fluctuations of the detection laser or a change 
in the bias magnetic field alter the photon distributions of both dark and bright states from 
shot-to-shot. If the state detection method depends on either form or average value of the 
photon distribution, the result of the state assignment will be biased by a change in the 
scattering rate. 



3.2. Threshold Detection Technique 



In this method, state discrimination is achieved by imposing a threshold o and assigning all 
detection events which yield a photon number higher than the chosen threshold to the bright 
state and all events with equal or lower photon numbers to the dark state. The advantage of 
this method lies in its simplicity and the fact that a single detection cycle suffices to determine 
the state the ion was projected onto. 

The state detection error for the bright state is given by 

r(<r+i,f B ) 



Pl = p(^<a) 



ifc=0 



For the dark state, it reads in the absence of depumping 



pI = p&>°) 



k=G+l 



kl Ol 

where r(<7, B, ) and y(<7, £ ) are the incomplete upper and lower gamma functions 



r(<x,l: 



dte-'t*- 1 



dte-'t a - 1 



(5) 



(6) 



(7) 



Since depumping is neglected the photon counts of the dark states correspond to the 
background, i.e. £,d = The depumping changes only the dark state and is again included 
by weighting the probability distributions with an exponential decay w(t) again 

-T _ poo _ 

dt w(t) p(k\$t) + dt w{t) p(k\&o) . (8) 



Pl=p(^>a)= £ 

k=a+\ 

The first term describes the case in which a depumping event took place during the detection 
time T. This effectively changes the mean value of the Poissonian photon distribution to 
4 = Roat + 7? b (t — t) ll27ll . The second term accounts for the case where only background 
counts are measured, i.e. the depumping took place after the detection. The expression can be 
further simplified to 

" a -(*+l) 

[y(fc+l,<f B -a)-y(fc+l,<^-a)] + 



J3 
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k=a+l 
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(9) 
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with the definitions /3 = (£;„ — t,^) ■ T/x and a = 1 — i. In Figure [3](a), the average detection 
error 

e th = (10) 

for the threshold technique is shown as a function of the detection time x and the threshold o. 
An increase in the detection time leads to smaller errors for the bright state, but the probability 
of a depumping event increases, resulting in larger errors for the dark state detection. On the 
other hand, a higher threshold decreases the dark state error, but larger parts of the bright 
state photon distribution are wrongly assigned to the dark state. For this reason, there is an 
optimum detection time for each chosen threshold which minimizes the threshold detection 
error. 

Among the main error sources of the threshold technique is the fact that the bright and 
dark state photon distributions overlap. We find that replacing the chosen threshold by an 
interval (o — 5, o + 8) which excludes all detection events in the overlapping region does 
not solve this issue. Instead, in this scenario, the distributions are truncated asymmetrically, 
especially with an average low number of detected photons, leading to a bias in the assigned 
state of the ion (see also Section [3T4b . 




Detection Time % (us) Detection Time x (\is) 

Figure 3. Theoretical detection errors in percentage (%). The scattering rates for both 
states are R B ~ 146.3 • 10 3 s _1 (R^ ~ 2.9 • 10 3 s _1 ) in this calculation. Part (a) depicts the 
expected error (Eq. [Tol l for the threshold technique, whereas part (b) shows the expected error 
for the ^-detection (black contour lines) as a function of the detection time and the threshold 
(see Section [3~4l for details). Additionally shown in part (b) is the approximate amount of 
the remaining statistics (%) according to Eq. \A.\\ (red contour lines). This value reflects the 
fraction of measurements that are not discarded by the post-selective filtering. 



3.3. Distribution-Fit Detection 

In this detection method, the full distribution of the photon histograms of both bright and 
dark states are taken into account for state discrimination. Prior to the experiments, two 
calibration photon histograms are taken: one where the ion is prepared in the | \) state by 
means of optical pumping (p c (fc|| B )) and a second one where the resonant Doppler-cooling 
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laser is either switched off or the ion is released from the trap and only background counts are 
observed (p c (fc||oo)) [20]. Here, p c (k\E,) corresponds to the calibrated probability distribution 
of measuring k photons given an average number of photons of f . The bright state amplitude 
a is then determined by a least-square fit to a superposition of both distributions 

p(k\a)=a-p<(k\$ B ) + (\-a)-pc(k\^) . (11) 

The expected detection error of this distribution-fit detection can be determined by Monte- 
Carlo simulations and is shown in Figure |4] as a function of the number of total detection 
events. As postulated by the law of large numbers, the error decreases with the square- 
root of the number of experiments. For our experimental parameters of typical 300-1000 
experiments, a value of 2-4% is expected, which is consistent with the standard deviation of 
the measurement result using the distribution-fit technique. 




'o 100 200 300 400 500 600 700 800 900 1000 
No. of Experiments 



Figure 4. Simulated detection error of the distribution-fit technique. The amount of scattered 
photons is determined in a Monte-Carlo type simulation assuming the ion is in an equal 
superposition of bright and dark states. The statistical fitting error of the least square fit of 
the photon histogram to Eq. [TTJ is plotted as a function of the number of experiments. The 
scattering rates are the same as in FigureO 

While the predicted error under the assumption of our typical experimental parameters 
is smaller than the error of the threshold detection, a sufficient number of measurements is 
required to fit the photon distributions. This detection method does not provide single-shot 
results. 

3.4. TZ-Detection 

In order to overcome the described limitations of the previous methods, two detection events 
which use the threshold technique are combined with an intermediate spin-flip and anti- 
correlations of both detection results are analyzed. The detection pulse sequence is shown in 
Figure 0(a). If both events yield the same state result, the detection is identified as false and 
discarded. For instance, assume the ion is in the bright state before detection, but is assigned 
to the dark state since the number of detected photons is below the threshold. The spin-flip 
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brings the ion into the dark state where it is detected for a second time. If this measurement 
(correctly) yields a dark state again, the whole detection process is discarded, since a false 
detection was recognized. Thus, the K-detection method effectively acts as a statistical filter 
for the measured data. 

The detection error for this scheme is determined by combining the errors of two 
threshold detection events and a finite fidelity of the intermediate spin-flip of 1 — £ lf , while 
neglecting discarded events. All possible paths of the whole detection process are shown in 
Figure [6l assuming that the ion is initially in the dark state. The decision tree for the bright 
state is similar. Summing over all possible events yields for the dark state error including 
depumping effects 

iC= (12) 



J'dt w(t) p{k > {p(^ <a) + e,(p(^ < a) -pfa < a))} 

+ dt w{t) p(k > a||oo)J {p(& < a) -e a (p(£ E < ct) < o))} 



J8 



— [y(*+l,&-a)-y(*+l,$»-a 



k=o+l 

■ {p(& < o)+e rf (p^ B < o)-p{&, < a))} 
+ e -*/T. r(°+J>& {p{ ^ < a)-e ri (p(^ < a) < a))} . 
On the other hand, if the ion is initially in the bright state, the error is given by 

Ke„ = < o) ■ {p(^ D > a)+£ rr (p(^ > o) -pfa >a))} . (13) 

The advantage of this method is that two small probabilities are multiplied, yielding an overall 
smaller error than the pure threshold detection method. In Figure [3](b), the average detection 
error of bright and dark states as a function of the chosen threshold o and the detection time 
T is shown. While the error is almost an order of magnitude smaller than for the threshold 
technique and tends to zero for very small detection times, it should be noted that this is at 
the expense of the usable data, also shown in the figure. For this calculation, the intermediate 
radio-frequency transition was assumed to have a fidelity of 98%, which is a conservative 
lower bound of the actual fidelity G8l . If necessary, the implementation of different, more 
robust techniques that are independent of the area of the radio-frequency pulse, such as rapid 
adiabatic fast passage [|29l or STIRAP 11301 can improve the fidelity even further. 

While the post-selective analysis for the ^-detection method reduces the detection error 
by discarding all photon counting events in the overlapping region of the histograms, events 
which undergo depumping during detection introduce a remaining bias in the resulting state 
amplitude. This type of error is already taken into account in the above calculation, but shall be 
analysed here in more detail. As described in Section l3T2l an asymmetric discarding of events 
with respect to both the bright and the dark states will introduce a bias towards the state with 
less discarded events. This bias is thus readily determined by comparing the difference in the 
probabilities of correctly detecting either the bright or the dark state. Neglecting the infidelity 
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Figure 5. (a) ft-detection scheme. Two identical threshold detection events are combined 
with an interleaved radio-frequency induced spin-flip. The detection result is only accepted 
if both detection events yield opposite results (green boxes). In the case of identical results, 
the event is discarded (yellow boxes), (b) 7T-Bayesian detection. Similar to the ft-detection, 
two Bayesian detection events are connected with an intermediate spin-flip. Here, a single 
Bayesian detection event consists of a W consecutive sub-bins of length t s . The detection 
result of each pulse is determined by the maximum likelihood formula (see e.g. Eq.fTol). 



of the interleaved spin-flip, it reads 

K = p{&><o)-p(&>o) . (14) 
In the case of the ^-detection, a sum over all possible correct detection events yields for the 



bias 



K = p($b > o) ■ (J^ dt w(t) p(k < a||oc) j 
+ p(^ D >o)-^ dtw{t)p(k<o\&)J 



-p{^<o)- P {^>o) . (15) 

Both biases b th and b K are plotted in Figure [7] as a function of the decay time T . Evidently, the 
threshold method introduces a systematic bias to the measurement result. The bias is shifted 
towards the dark state with increasing decay time and levels at approximately 20 % for T — > 00, 
for the parameters described previously. In contrast to that, the ^-detection method operates at 
a smaller bias and culminates in a bias-free detection for systems with large decay times. This 
is a result of the symmetric approach of comparing two consecutive and inverted detection 
events. The decay time of the 25 Mg + qubit system being approximately T w 60 /is, the %- 
detection method is estimated to have a bias of —5.2%, whereas the bias of the conventional 
threshold method is estimated to be 13.1 %. 

3.5. Bayesian and K-Bayesian Detections 



In the previously discussed methods, the collected number of photons integrated over the 
detection time T was considered. This neglects all information which could be gained from 
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Figure 6. Decision tree for ^-detection. Assuming the ion to be initially in the dark state, 
all possible paths of both detection events with their corresponding probabilities are shown. 
Also included is a possible error in the spin-flip by assuming a finite fidelity of 1 — Eg for a 
successful inversion. The first branches are marked for clarity: I) Off -resonant excitation II) 
Threshold detection III) Spin-flip. Circles correspond to the actual state of the ion, whereas 
rectangular boxes denote the detection result (1 for the bright and for the dark state). The 
color coding distinguishes between correct detection events (green), discarded events (yellow) 
and detection errors (red). 



the arrival times of the photons. This information can be included by dividing the detection 
interval into N sub-intervals of length t s . The state of the ion is then inferred from the time- 
resolved series of detected photons where i = 1, . . . , N. The likelihoods p({£jj}|B,D) of 
observing a particular series of photons given that the ion is in either bright (B) or dark states 
(D) are calculated and the ion is assigned to the state with the highest likelihood lfT3ll2TTl . 
The likelihood of the photon series originating from the bright state is given by 



N 



p({&>iB) = rb(&iso 



k=0 



In case of the dark state, it reads 



-h/T 



~Nt s /T 
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(16) 
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Figure 7. Calculated bias b lh 7[ of the detection methods as a function of the decay time T. The 
black dotted line is a guide to the eye and corresponds to the zero bias line. The red dashed 
line reflects the case of the 25 Mg + system with a decay time of approximately T ~ 60/Xs. The 
calculation was done for a detection time of 10/Xs and a threshold of a = 0. The scattering 
rates R B D are the same as in Figure [3] 



rf», = 1 - PCD|{&}) =PCB|{6» = ., ft ,^ , (18) 



where it was assumed that depumping events only happen at the end of a sub-interval. The 
detection error is determined by Bayes' theorem, i.e. the probability that the bright state is 
wrongly deduced to be dark reads lfP3l 

p(mp)+p({tt}\B) 

and vice versa for the dark state. This can be used to determine the detection error in real time 
and stop the measurement as soon as the desired error level is reached EE]. 

In an identical way as for the threshold detection method, the Bayesian detection can 
be extended to a ^-Bayesian detection, as shown in Figure |5](b). Here, two detection events 
which use the maximum likelihood to determine the state of the ion are combined with an 
intermediate spin-flip. This again serves as a statistical filter and discards correlated events. 
The detection error is similarly calculated. For instance, while neglecting the imperfect 
fidelity of the spin-flip, the error of the bright state detection reads 

pL y ,J{^ ] ^ 2) }) =P(V\{& 1) }) -P(B|{^ 2) }) • (19) 

In Figure [8](a), calculations of the average detection error e th „ of both the threshold 
and ^-detection are shown in comparison as a function of the detection time, illustrating the 
effect of depumping on the detection error. While there is an optimal point for the threshold 
technique due to an increase in error from depumping for longer detection times, the Tl- 
detection yields smaller single-shot errors at the expense of the number of used events, as 
shown by the dashed blue curve. Consequently, the statistical error which results from the 
quantum projection noise increases as the number of used events decreases OTl . This is 
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shown by the plotting the 95% confidence interval £ th <K ± 1.96a/ a(l —a)/n (shaded areas in 
the plots), where a represents the bright state amplitude, similar to Eq. [TTJ The number of 
experiments is n = 1000 for the threshold technique and is scaled by the percentage of used 
events for the ^-techniques. Furthermore, we use a = 0.5 as a conservative estimate, owing 
to the maximal quantum projection noise if the ion is in an equal superposition of bright and 
dark state. 

Although the ^-detection discards many events, the resulting statistical spread still yields 
a smaller overall detection error than the threshold technique. Nevertheless, it should be noted 
that applying the ^-detection in a quantum algorithm with several consecutive detection steps, 
the algorithmic sequence needs to be repeated in case of a single inconclusive result. It is 
worthwhile mentioning that the non-zero single- shot error of the threshold technique leads to 
a bias that fundamentally does not average to zero, whereas the statistical uncertainty can be 
improved by increasing the number of measurements. 

The same holds for the Bayesian techniques. In Figure[8](b), Monte-Carlo simulations of 
both Bayesian and ^-Bayesian technique are shown in comparison. The threshold technique 
yields similar results to the Bayesian method for short detection times. The observation of the 
photon series only improves for detection times longer than 20 /is. Applying the ^-method 
to the Bayesian scheme additionally decreases the detection error for short detection times at 
the expense of the used events. 

Choosing an optimal threshold for each detection time for both threshold and ^-detection 
decreases the error for both methods further. Although smaller single-shot errors are achieved 
for longer detection times with an optimal threshold, there is a trade off by discarding more 
events with the ^-technique. 
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Figure 8. Theoretical error for different detection methods as a function of the detection 
time. The parameters are the same as in Figure [3] For the Bayesian method, a Monte-Carlo 
simulation with 2000 events at each point was used to determine the error. The dashed lines 
correspond to the amount of remaining statistics for the 7T-type methods. A threshold a = is 
chosen. The shaded areas reflect the expected additional statistical error for 1000 experiments, 
(a) Threshold and ^-detection methods, (b) Threshold, Bayesian and 7T-Bayesian methods. 
For the simulation of the Bayesian techniques, a sub-bin time of 2/is was used. 
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4. Robustness and Sensitivity of Detection Schemes 

A major difference among the detection methods is their sensitivity to the chosen threshold 
and to fluctuations in the atomic scattering rate from technical noise. While the Bayesian and 
distribution-fit detection methods are inherently threshold free, the threshold and ^-detection 
schemes may be affected by choosing the wrong threshold. This is discussed in detail in 
Appendix A 

Fluctuations of experimental parameters, such as the power and frequency of the 
employed detection laser or the strength of the external magnetic field, are a major 
experimental challenge since they change the fluorescence rate of the atom and introduce 
further errors into the measurement process. Thus, it is highly desirable to aim for robust 
detection methods. Here, sensitivity with respect to power fluctuations was measured for the 
different detection techniques. 

In order to simulate power fluctuations, the ion was initialized in an equal superposition 
of the 1 1} and | f) states by means of a n/2 radio-frequency pulse. The power of the detection 
laser was scanned over ±5 dB around its optimal value of dB (corresponding to saturation 
intensity), as shown in Figure [9l At each point, 3 x 250 measurements were taken. Since 
a change in laser power alters the photon distributions, one expects the detection methods 
that depend on the threshold to be especially sensitive towards power changes. This is 
confirmed in the experiment. While the threshold, distribution-fit and Bayesian techniques 
tend to over- and underestimate the actual 1 1) population of the ion by 30%, both ^-detection 
techniques show no significant dependency on power fluctuations, since wrongly assigned 
events are systematically discarded. Given a typical detection time of 10/is, the Bayesian 
and the threshold techniques yield comparable results, which is consistent with the fact 
that their individually expected error is similar, as shown in Figure [8](b). Additionally, the 
expected statistical error due to the finite number of measurements, similarly to Figure [8l is 
indicated by the shaded areas. A linear fit to the interval between -1 dB and +1 dB yields 
the following changes of the amplitude for the different techniques: (a) threshold technique: 
0.032 ± 0.01 1/dB, (b) ^-detection: -0.009 ± 0.016/dB and (c) distribution-fit technique: 
0.037 ±0.015/dB. 

It is important to stress that the systematic bias found in the conventional detection 
methods dominates the overall error in comparison to the statistical error. Although the 
statistical error can be decreased by obtaining more measurements, the systematic bias is 
only corrected in the ^-detection methods. 

5. Conclusion 

We demonstrated a novel qubit detection technique which combines two detection events 
with an intermediate state inversion. The expected anti-correlation of both detection events is 
used in a post-selective statistical filtering. As a consequence, a higher detection fidelity 
is achieved at the expense of the number of usable events. Although this results in an 
increase of the statistical error due to quantum projection noise, the systematic bias error 
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Figure 9. Robustness with respect to power fluctuations for all detection methods. The ion 
is initialized in an equal superposition of bright and dark states and the nominal detection 
power is dB, corresponding to the saturation intensity. Detection methods that depend on the 
threshold experience large systematic deviations from the initial state result. Both ^-detection 
techniques are almost independent of the laser power at the expense of usable events. Shaded 
areas depict the statistical error as in Figure [8] The detection time for the threshold detection 
and both threshold events of the ft-detection was T = 10/is. In the case of the Bayesian 
detections five sub-bins with a length t s —2jj,s were used. 



of the conventional threshold technique is in comparison significantly larger and does not 
average to zero. We compared the ^-detection to other methods in terms of their sensitivity 
to fluctuations of experimental parameters. In particular, the ^-detection method shows no 
significant dependence on the detection time, the chosen threshold or power fluctuations of 
the detection laser, which makes it a rather robust experimental tool. 

Furthermore, especially in a regime of low photon count rate, the overlap of the photon 
histograms for bright and dark state imposes a limitation on the achievable contrast of Rabi 
oscillations to the conventional threshold technique, whereas almost full contrast is observed 
using the ^-detection (see Appendix A). 

Another advantage of this method is its simplicity, requiring only an rf-induced state 
inversion for qubits encoded in hyperfine ground states. This technique applied to optically 
encoded qubits would require a narrow-bandwidth laser for coherent manipulation, which is 
often available in ion trap quantum computing and precision spectroscopy experiments. 

The ^-detection method could also be applied to multiple qubit systems where the qubits 
are individually accessible for detection, as e.g. in ion traps where multiple ions are spatially 
resolved ll32ll . Furthermore, the concept of this detection scheme is very general and can 
be applied to a multitude of systems for which non-destructive detection methods can be 
implemented, such as neutral atoms in deep optical lattices ll33l l34l . measurements of the 
state-dependent reflection and transmission signals of molecules or atoms in cavities [|22ll23l . 
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and nitrogen vacancy centers 1(35113611 . 
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Appendix A. Threshold Sensitivity 

We investigated the sensitivity of the detection methods to the threshold by observing the 
contrast of radio-frequency induced Rabi oscillations. Figure [ATI shows a single experimental 
data set analyzed with different detection techniques. In subplot (a), Rabi oscillations are 
shown for the threshold and the distribution-fit method for a threshold of 1 photon. Clearly, 
the contrast resulting from the threshold technique is degraded, owing to the non-optimal 
threshold; subplot (b) depicts the analysis of the data with the ^-detection compared to the 
threshold method for the same threshold of 1 photon. The ^-detection yields almost unity 
contrast by discarding over 50% of the events. These correspond to the events that were 
wrongly assigned to the dark state by the threshold method. This effect can be understood 
by considering the amount of statistics used by the ^-detection. Neglecting the fidelity of the 
radio-frequency and the depumping effect for simplicity, the percentage of used statistics is 
given by 

P^ at = l-p(^<a)-p(^ 00 > a)+2-p(^<a)-p(^> a) . (A.l) 

In a regime where the threshold is so high that the error of the dark state detection is negligible, 
the remaining statistics correspond to p n stat « 1 — p(£„ < <j) . On the other hand, applying the 
threshold technique in this regime, the amplitude of the Rabi oscillations is determined by the 
error of the bright detection, i.e. A th = p n sat which is confirmed by the measurements in Figure 

ED(b). 

Similarly, the 7T-Bayesian detection filters correlated events and yields an almost unity 
contrast as opposed to the pure Bayesian detection. Both are depicted in Figure lATI Cc). In 
subplot (d), the dependence of all methods is summarized by plotting the fitted contrast 
of the Rabi oscillations as a function of the threshold. While the distribution and both 
Bayesian techniques are fundamentally threshold free, the threshold method experiences a 
strong decrease of the contrast with increasing threshold. On the other hand, the ^-detection 
discards many events and the resulting contrast is almost independent of the threshold. 
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